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ABSTRACT 

We report the discovery of a new ultra-faint globular cluster in the constella- 
tion of Ursa Minor, based on stellar photometry from the MegaCam imager at the 
Canada- France-Hawaii Telescope (CFHT). We find that this cluster, Munoz 1, 
is located at a distance of 45 ± 5kpc and at a projected distance of only 45' 
from the center of the Ursa Minor dSph galaxy Using a Maximum Likelihood 
technique we measure a half-light radius of 0'.5, or equivalently 7pc and an el- 
lipticity consistent with being zero. We estimate its absolute magnitude to be 
My = —0.4 ± 0.9, which corresponds to Ly = 12Oi 65 °L and we measure a 
heliocentric radial velocity of —137 ± 4 km s -1 based on Keck/DEIMOS spec- 
troscopy. This new satellite is separate from Ursa Minor by ~ 30 kpc and 110 km 
s _1 suggesting the cluster is not obviously associated with the dSph, despite the 
very close angular separation. Based on its photometric properties and structural 
parameters we conclude that Munoz 1 is a new ultra-faint stellar cluster. Along 
with Segue 3 this is one of the faintest stellar clusters known to date. 
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Introduction 



(SDSS; 


York et al. 


2000), seventeen new ' 


Willman et al. 


12005: 


Zucker et al. 


2006a 


b 



The majority of these objects correspond to a new class of ultra low-luminosity dwarf galax- 



ies, based on their kinemat ic and metallicity properties (e.g. JMunoz et al.ll2006l : Martin et al. 
20071 : ISimon fc Gehall2007l ) , with luminosities ranging from My ~ — 8 for the br ighter of these 
objects down to an extreme My ~ —1.5 for the faintest (jMartin et al.l 1200a ). Given their 



extremely low stellar contents, they are commonly referred to as the Ultra-Faint Dwarfs 
(UFDs). A few other systems are clearly low l uminosity globular c lusters, which include 
Koposov 1 and 2 (My = —1 and —2 respectively; IKoposov et al.ll2OO70 and the extreme case 



of Segue 3 (My ~ 0.0; iBelokurov et al 



20101 : iFadelv et al.ll201lh 



Despite the enormous success of satellites searches in SDSS, this survey covers only 
about a fourth of the sky down to a magnitude limit of r ~ 22.5. In practice this means 
that the faintest objects can only be detected out to distances smaller than a few tens of 
kiloparsecs (see Fig. 9 in iTollerud et al.l 120081 ). and thus the majority of the Milky Way's 
virial volume remains unsearched for UFDs and other faint stellar systems. 

In this Letter we present the discovery of a new extremely low-luminosity, outer halo 
stellar cluster, in the outskirts of the Ursa Minor dwarf spheroidal (dSph) galaxy. This 
discovery was made just outside the SDSS footprint using significantly deeper photometry 
which allowed us to probe farther into the Milky Way halo. In §2 we detail the data-set and 
the discovery. In §3 we present the results from the photometric analysis and in §4 we show 
radial velocity data recently obtained. Finally, in §5 we discuss and summarize our results. 



2. Data and Discovery 

The discovery of this new system was made serendipitously while analyzing photometric 
data for the Ursa Minor dSph galaxy taken with the MegaCam imager at the Canada- France- 
Hawaii Telescope (CFHT). MegaCam is a wide-field imager consisting of 36 2048 x 4612 pixel 
CCDs, covering almost a full lxl deg 2 field of view with a pixel scale of 0" .187 pixel -1 . These 
data were taken as part of a larger program aimed at obtaining deep wide-field imaging of 
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all bound stellar over-densities in the Milky Way halo beyond 25kpc (R. R. Munoz et al. 
2012, in preparation). Ursa Minor was observed on the nights of UT June 12-14 and July 
07-13, 2010 under dark conditions with typical seeing of 0.7 — 0".9. Four different, slightly 
overlapping fields were observed for a total area coverage of nearly 2 x 2deg 2 . In each field, 
the center of the dSph was placed in one of the corners so that, when combined, the galaxy is 
located at the center of the covered area. We obtained six dithered exposures of 360 seconds 
in both the g— and r— band. A standard dithering pattern was used to cover both small and 
large gaps present between the chips. 

MegaCam i mages are delivered to the us er pre-processed by the CFHT team using the 
"Elixir" package (jMagnier &: Cuillandrdl2004l ). Subsequent photometric measurements were 
carried out using first DAOPHOT/Allstar, and later runni ng the ALLFR AME package on 
the processed frames, following the procedure outlined in Istetsonl Jl994h . Finally, astro- 
metric solutions were calculated using the freely available SCAMF0 code, and photometric 
calibra tion was carried out by direct comparison w ith data from the SDSS-Data Release 7 
fDR7; lAbazaiian fc Sloan Digital Sky Survey! 120091 ). For more details on the method see 



Munoz et all (120101 ) . 



As seen in Figure [T] and the left panel of Figure [21 we discovered a centrally concentrated 
over-density of stars 45' to the south-west of Ursa Minor. Based on the structural parameters 
determined below, we presume this object to be an ultra-faint globular cluster and therefore 
name it Munoz 1. In §5 we discuss in more depth our reasoning for this classification. 
Figure [2] also shows the color-magnitude diagram (CMD) of the central region of Munoz 1 
where a distinct main-sequence turn-off is observed. Assuming the stars in the red-giant 
branch (RGB) region are members of the cluster we fit an isochrone at a distance of 45 kpc, 
an age of 12.5 Gyr and [Fe/H]= —1.5. We note that this fit is very tentative given the very 
low number of stars that belong to the system. The right panel of the figure compares the 
CMDs of UMi and the new object showing a clear difference in distance modulus. 



3. Structural Properties 



We carry out a ma ximum-likelihood an alysis of the phot ometric data for M unoz 1 
following the method of iMartin et al.l ( 120081 ) as described in iMunoz et al.l (|2010j). This 
method starts by assuming a shape for the underlying light distribu tion; in this case we have 
tried three of t he most common ly used profiles for UFDs, a King (jKing 1962 ). exponential 
and Plummer ( jPlummerl Il91ll ) density laws. We then fit simultaneously the structural 



See http: /www. astromatic.net/software/scamp. 
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parameters of the object, i.e., scale length (half-light radius or King core and tidal radii 
depending on the profile), coordinates of the center of the cluster, ellipticity, position angle 
and background density. To estimate parameter uncertainties we carry out a bootstrap 
analysis of 10000 realizations of the photometric data. 

Table 1 shows the resulting parameters. We obtain a half-light radius of 7.1 ±2.1 pc using 
both an exponential and Plummer profile, whereas a King profile yields r core = 4.5 ± 1.9 pc 
and r tidal = 65 ± 28 pc. 

Figure |3] shows the number density profile for Munoz 1. We have overp lotted the three 
best-fit profiles, all of which give a reasonable description of the light distribution. The right 
panel of this figure shows a density contour map of the object, where the contours represent 
2, 3, 5, 8, 12 and 18 — a levels above the background density. These contours show no 
evident elongation or tidal features. A pplying a bootstrap analysis where we generate 10000 
realizations of the photometric data ( iWalsh et al.l 120081 ; iMunoz et al.ll2010l ). we determine 
that the asymmetries observed in the outer parts, toward the south and west directions, are 
not statistically significant due to the very low number of stars. 

In addition to the structural parameters, we estimate t he ab solute magnitude of the 
new object. We follow the method described in lMunoz et al.l ( 120101 ) which relies in the total 
number of stars tha t belong to the cluster and not on the sum of their fluxes. As shown by 
Martin et al.l ( 120081 ) the low number of stars identified in the ultra-faint systems, especially 
true in the case of Munoz 1, makes traditional methods of adding the individual fluxes of the 
member stars very sensitive to the inclusion (or exclusion) of potential members, especially 
at brighter magnitudes. To alleviate shot noise issues we use an alternative method based 
on a model stellar population. For this we use a theoretical luminosity function that best 
describes the p hotometric propert ies of the cluster, in this case a 12.5 Gyr old population with 
[Fe/H]= —1.5 ( iDotter et al.l 120081 ) . We then integrate the LF to obtain the total flux down 
to a given magnitude limit. The last step is to scale this flux using the total number of stars 
that belong to the cluster down to the same magnitude limit. We estimate the corresponding 
uncertainty through a bootstra p analysis . Thi s yields My = —0.4 ± 0.9 mag assuming a 



Chabrier initial mass function (jChabrierl l200ll ). which corresponds to Ly = 120tg5°L Q . 



jy — -l^w_ 65 

Only one old star cluster has a published luminosity lower than Munoz 1, the Segue 3 stellar 
cluster with a total luminosity of My = —0.0 ± 0.8 mag (Fadely et al. 2011), although given 
the large uncertainties in both measurements it is unclear which cluster is actually fainter. 
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4. Spectroscopic Properties 



Spectroscopic data of photometrically selected Munoz 1 candidate stars were taken 
on UT May 28, 2011 with the Keck II 10-m telescope and the DEIMOS spectrograph 
( Faber et al. I l2003h . One Keck/DEIMOS multi-slit mask was observed with the 1200 line mm 1 grating 
covering a wavelength region 6400 — 9100A with a spectral dispersion of 0.33A. The mask 
was observed for 7800 seconds under relatively poor seeing conditions which varied between 
1.5 — 2". Spectra were reduced using a modified version of the spec2d software pipeline 
(version 1.1.4) developed by the DEEP2 t eam, and radial velocities were determined using 
th e software d e scribe d by ISimon fc Gehal (120071 ). For additional details, we refer the reader 



to 



Geha et all (120091 ) . 



Radial velocities were successfully measured for 24 of the 47 extracted spectra. Stars for 
which we could not measure a velocity were primarily very low S/N spectra of faint objects 
near the expected main sequence turn-off of Munoz 1 (r ~ 22 — 22.5). 

The resulting velocity distribution is shown in Figure HI along with the corresponding 
spatial and color-magnitude distributions. While the velocity histogram does not show a 
clear significant peak that could be interpreted as the velocity signature of Munoz 1, if we 
restrict ourselves to the most centrally concentrated stars, what appears as a cold peak is 
indeed discerned. Using the four stars closest to the derived center of the cluster, plus a 
fifth star possibly associated with it based on its radial velocity, distance to the cluster and 
position in the CMD, we estimate a mean velocity of —137 ± 4 km s _1 . Given the small 
number of stars and large velocity uncertainties, we estimate only a one sigma upper limit 
on the velocity dispersion of o < 4.7 km s -1 . Given the size and lumin osity of Munoz 1 , 
we predict its velocity dispersion to be between 0.2 — 0.3 km s" 1 based on IWolf et al.l (120101 ) 
and a mass-to-light ratio consistent with old stars only. Our measured velocity dispersion is 
consistent with this value, but we cannot rule out a scenario where the mass of Munoz 1 is 
significantly larger than its stellar mass. 

Two other stars with radial velocities in the vicinity of this peak do not live near the 
CMD of Munoz 1, nor are they close to the center of the object, and therefore we consider 
them less likely to be associated with the cluster. We also observe three stars consistent 



with the published radial velocity of Ursa Minor of —245 km s _1 (IMunoz et al.l 120051 ). These 
stars are less spatially concentrated than the Munoz 1 stars and yield a mean velocity of 
—253 ± 6 km s -1 , in good agreement with the observed value for the dSph galaxy. 

A subset of stars in our Keck/DEIMOS kinematic sample have sufficient S/N (at least 
10 A ) to determine th eir metallicit i es. W e determine [Fe/H] based on fitting synthetic 
spectra, as described in iKirby et al.l (120081 ). This medium- resolution technique allows us 
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to calculate T e g- and [Fe/H] by matching an observed Keck/DEIMOS spectrum to a grid of 
models with different T e g and [Fe/H]. This is in contrast to most high resolution abundances, 
which are based on equivalent widths and require resolved, unblended lines. Our spectral 
resolution is not sufficient for such type of analysis. We note, however, that our well- 
characterized noise array enables us to calculate realistic error bars. 

We measure metallicities for two stars kinematically associated with Ursa Minor and 
one star associated with Munoz 1. The two Ursa Minor stars have metallicities [Fe/H] = 
-2.78±0.60 and -2.30±0.68, while the Munoz 1 star has [Fe/H]= -1.46±0.32. These results 
are consistent with the isochrone-derived metallicity. Due to the paucity of Fe II absorption 
lines in cool RGB spectra, we are unable to calculate surface gravity from the spectroscopic 
data. Hence, to estimate [Fe/H] we had to fix log g from the best-fit isochrone. The code also 
outputs a value for T eff , which can be compared to the value obtained using the de- reddened 
CFHT photometry. For the case of the cluster star, we obtain a T c g = 5160 ± 320 K, in good 
agreement with the photometric value of T e g = 5290 K. 



5. Summary and Discussion 



We report the discovery of a new ultra-faint stellar system, Munoz 1, in the Ursa Minor 
constellation, located only 45' from the center of the Ursa Minor dSph. The discovery 
was made serendipitously while analyzing photometric data for Ursa Minor, taken with the 
wide-field MegaCam i mager on the CFH T. We visually fit a 12.5 Gyr isochrone from the 
Dartmouth database (IDotter et al.l 120081 ) . with a metallicity of [Fe/H]= —1.5 at a distance 
of ~ 45kpc. While the close angular distance between these two Galactic satellites might 
suggest an association, we deem it unlikely. We have estimated the escape velocity from 
Ursa Minor at the p rojected distance of Munoz 1, i.e., lkpc, assuming the dSph is a point 
mass of 1 x 1O 8 M ( jWolf et al.ll2010l ). We obtain t> esc ~ 30 km s -1 , much smaller than the 
~ 110 km s~ l of radial velocity difference. Moreover, at 30kpc, the distance separating the 
cluster from the dSph along the line of sight, t> esc decreases down to only a few km s _1 . 



Using a maximum-likelihood technique, following the method of iMartin et al.l ( 120081 ). 
we have simultaneously derived structural parameters along with background density for 
Munoz 1. Assuming different density profiles we have derived consistent half light radii of 
r h,ex P = 7.1 ± 2.8 pc, r h P = 7.1 ± 2.1 pc for exponential and Plummer profiles respectively, 
and a King core and tidal radii of r c = 4.5 ± 1.9 and r t = 65 ± 28 pc respectively. The 
total luminos i ty of Munoz 1, on the other hand, was derived using the method described in 
Munoz et al.l ( 120101 ). which relies in the total number of stars in the object instead of their 
individual fluxes, obtaining an absolute magnitude of My = —0.4 ± 0.9, equivalent to Ly = 
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120lg5°-^©- Based on its photometric properties, we argue that this new Galactic satellite 
corresponds to an extremely faint globular cluster. Since the discovery of the UFDs, there has 
been some debate in the literature about the classification of stellar over-densities based solely 
on photometric information. It has been argued, for instance, that objects like Segue 1 and 
Coma Berenices should be considered large globular clus ters based on the as sumption that 
galaxies cannot have scale lengths smaller than ~ 100 pc (IGilmore et al.ll2007l ). On the other 
hand, spectroscopic information, i.e., velocity dispersions and metallicities, have been used to 
argue that despite their sizes, these syst ems correspond to dark matter dominated galaxies 
(jSimon &: Gehal 120071 ; iGeha et al.l 120091 ). For Munoz 1, the spectroscopic information is 
insufficient to determine unambiguously whether the object is dark matter dominated or not. 
However, this satellite is significantly smaller than any of the known dwarf galaxies, including 
the UFDs, it lies within the size distribution of other Milky Way globula r clusters with 



struc tural parameters similar to other outer halo clusters like Palomar 13 ([Bradford et al. 
201ll ). and while in principle it could be the smallest known galaxy, without a reliable 



measurement of its velocity dispersion of metallicity spread, we view it as much more likely 
that it is a globular cluster. The luminosity derived for this object, makes it one of the 
faintest globular cluster discovered, o nly nominally surp assed by Segue 3, which has an 
absolute magnitude of My = 0.0 ± 0.8 (iFadely et al.ll201ll ). Given the large uncertainties of 
both measurements, it is unclear which one of them is indeed fainter. The density profile 
of this object is well fitted by typical cluster density profiles showing no obvious hints for 
significant tidal disruption. 

The discovery of this cluster was made in the context of a larger program aimed at 
obtaining deep and wide-field photometry for all outer halo satellites regardless of mor- 
phological classification. At the time of this publication, we have imaged roughly 40 deg 2 
down to a magnitude limit of r ~ 24.5 — 25 depending on the object. As pointed out by 
Tollerud et al.l (120081 ). surveys deeper than SDSS are bound to uncover fainter and more 
distant ultra-faint systems. Naively extrapolating the discovery of at least one ultra-faint 
cluster in our surveyed area to the entire sky, yields potentially a thousand similar objects 
still uncovered. While this number is not intended as a firm prediction given its large un- 
certainty, it foreshadows a much more complex picture of the globular cluster system in the 
outer halo than we have today. 
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Table 1. Parameters of Munoz 1 



Parameter 


Mean 


Uncertainty 


a 0:exp (h m s) 


15:01:48.02 


±13 


0Q,exp (d m Sj 


+oo:oo:U ( .6 


_L all 


^/i,ea;p (arcmin) 


U.4y 


in i n 

±u.iy 


?Vsp (pc) 


( .1 


1 o o 

±z.o 


r^ )P (arcmin) 


U.49 


±U.lo 


rh,p (pc) 


7.1 


1 o -i 

±2.1 


r c (arcmin) 


15" 


±9" 


r c (pc) 


4.5 


±1.9 


r t (arcmin) 


3.6 


±2.2 


r t (pc) 


65 


±28 


M v (Chabrier) 


-0.4 


±0.9 


Lv (L ) 


120 


±160,-65 


/ioy (mag arcsec -2 ) 


26.3 


±1.6,-2.1 


t> r (km s -1 ) 


-137 


±4 


d (kpc) 


45 


±5 




Fig. 1. — 1'.2 x 1'.2, r— band view of Munoz 1. The image is roughly centered on the cluster. 
North is up, east is to the left. 
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Fig. 2. — left: Star count map of the Ursa Minor dSph field. The primary object in this 
field is the dSph galaxy. Munoz 1 lies 45' away from the center of the field in the south-west 
direction, near (—1.8°,— 0.25°). middle: g vs g — r CMD of stars in the Munoz 1 region 
(within 2' of t he measured cente r). The best fit isochrone is overplotted, corresponding to 
a Dartmouth ( jDotter et al.l 120081 ) . 12.5 Gyr, [Fe/H]= —1. 5 at a distance of ~ 45kpc using a 
reddening value of E(g — r) = 0.027, adopted from the ISchlegel et al.l ( 119981 ) maps, right: 
CMD of Ursa Minor stars (small black points) overplotted on Munoz 1 stars (large red 
symbols). Red and blue solid lines represent best fit isochrones for both the Ursa Minor 
dSph and Munoz 1 respectively. 
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Fig. 3. — Left: Number density profile for Muiioz 1. Plummer, King and exponential 
profiles have been fitted to the data using a maximum likelihood method. The dot-dashed 
line represents the measured background density. Right: Iso-density contours for Munoz 1. 
Contours level correspond to 2,3,5,8,12 and 18— a level over the background density. 




Fig. 4. — Keck/DEIMOS spectroscopic coverage of Munoz 1, showing the spatial distribution 
(left), co lor- magnitude diagram (middle) and velocity distribution of targets (right). Red 
symbols in each panel indicate stars possibly associated with Munoz 1, while blue symbols 
are stars likely associated with Ursa Minor. 



